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ABSTRACT 

A bright, nearby edge-on starburst galaxy NGC 253 was studied using the Suzaku, XMM and 
Chandra X-ray observatories. We detected with Suzaku and XMM complex line structure of Fe K, 
which is resolved into three lines (Fe I at 6.4 keV, Fe XXV at 6.7 keV and Fe XXVI at 7.0 keV) around 
the center of NGC 253. Especially, the Fe I and Fe XXVI lines are the first clear detections, with a 
significance of >99.99 % and 99.89 % estimated by a Monte Carlo procedure. Imaging spectroscopy 
with Chandra revealed that the emission is distributed in ~60 arcsec 2 region around the nucleus, 
which suggests that the source is not only the buried AGN. The flux of highly ionized Fe lines can be 
explained by the accumulation of 10-1000 supernova remnants that are the result of high starforming 
activity, while the Fe I line flux is consistent with the fluorescent line emission expected with the 
molecular clouds in the region. 

Subject headings: galaxies: individual (NGC 253) — galaxies: ISM — galaxies: starburst — galaxies: 
star formation — X-rays: galaxies 



1. INTRODUCTION 

X-ray emission from the central region of a starburst 
galaxy is known to be complex because of blending of 
various kinds of point sources and diffuse emission. Al- 
though point source emission typically dominates the 
hard band, in some starburst galaxies including NGC 
253, the emission line from He-like Fe XXV is observed, 
sugg esting the existence of a high temperature plasma 
Ce.g.lPietsch et al.ll200lL iBoller et al1l200l iRanalli et al.1 
2008). Supernova remnants (SNRs), associated with 
starburst activity, are considered to be a candidate of 
the origin. Moreover, two of the g alaxies (NGC6240 
and M82) exhibit Fe I fluores cent line (jBoller et al.ll2003t 
Strickland fe Heckmanl[2007i ). The coexistence of highly 
ionized and neutral Fe lines suggests the coexistence of 
hot and cold gases. Thus, this complex Fe K line struc- 
ture is informative on the origin of the emission and the 
gas structure of starburst galaxies. The origin of the 
neutral line is, however, still poorly known. 

NGC 253 is a nearby bright edge-on starburst galaxy 
with the star for mation rate of 1.4-9.5 M Q yr _1 
(jGrimm et al.ll2003() . Its large apparent diameter makes 
it a suitable target to study the spacial distribution of 
emission features. Chandra suggested that the X-ray 
emission in the nuclear region is composed of a low- 
temperature thermal plasma and a heavily absorbed hard 
component that is expressed by a power law plus emis- 
sion from a photoionized plasma emitting Fe IV-Fe XV 
and Fe XXV (jWeaver et al.ll2002D . The hard component 
is int erpreted as a buried AGN (|Muller- Sanchez et al.l 
I201CO . The spectrum around the nuclear region could be 
modeled with three-temperature thin thermal plasmas 
including a high (^6 ke V) temperature plas ma responsi- 
ble for the Fe XXV line (|Pietsch et al.ll2001l ). 

In this Letter, we report on the study of the nuclear re- 
gion of NGC 253 using Suzaku and Chandra, with a par- 
ticular focus on the Fe line complex. Suzaku has a good 
energy resolution to resolve the Fe line complex with a 
low background, while a sub-arcsecond spatial resolution 
of Chandra enables us to argue spatial distribution. We 



also show results by XMM-Newton as a complementary 
analysis. Throughout this paper, we adopt the d istance 
to NGC 253 of 3.4 Mpc (palcanton et al.l [20091 ) corre- 
sponding to 16 pc arcsec" 1 . We assum e the s olar abun- 
dance tabulated by [Anders^Grewiii dHH). HEAsoft 
version 6.11, SAS 10.0.0, CIAO 4.3 and XSPEC 12.7.0 
were used for data reduction and spectral analysis. For 
all fitting processes, we adopted C statistics and binned 
spectra so that each bin contains > 10 counts after a 
background subtraction. We quote an error as a 90 % 
significance level of a single parameter in all tables, de- 
termined by Monte Carlo simulations and most of errors 
are consistent with those obtained by AC-statistic of 2.7. 

2. OBSERVATIONS & DATA REDUCTION 

NGC 253 was observed by Suzaku (seq:805018010) 
with an exposure time of 101 ks. We analyzed only XIS- 
FI (XIS0 and XIS3), which has a lower background than 
XIS-BI (XIS1) at > 5 keV to avoid a systematic uncer- 
tainty associated with a simultaneous fitting of different 
sensors. The redistribution matrix files (RMFs) and an- 
cillary response files (ARFs) were created with th e xisrm- 
fgen and xissimarfgen ftools (Is hisaki et al .1120071) , respec- 
tively. Since the bright region is within several tens of 
arcseconds, which is much smaller than the Suzaku point 
spread function (PSF; 1'.8 as a half-power diameter), we 
assume the source to be a point-source when creating the 
arf files. The accumulation of dark Earth observations 
was used to estimate the non X- ray background (N XB) 
spectra, via the xisnxbgen ftool (jTawa et al.ll2008|) . We 
subtracted NXB, the dominant background above 5 keV, 
but not the X-ray background, which was evaluated by 
an offset observation next to our field as to be less than 
3 % of the flux in R0. The spectra and responses of 
XIS0 and XIS3 were merged using the ftools mathpha, 
marfrmf and addrmf. 

We utilized PN data from two XMM observations (Ob- 
sID:0125960101 and 0152020101) with the effective expo- 
sure time of 31 and 54 ksec after removing flared time in- 
tervals which is over 0.35 counts/sec above 10 keV. RMS 
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Fig. 1. — Images taken by Suzaku XIS3 (left), XMMPN (middle) and Chandra ACIS (right) in 4.0-10 keV. The red ellipse and magenta 
cross mark correspond to an optical disk and the center of the Chanda image. Green circles and ellipses are regions used for extracting 
spectra. 



and ARF files were made by the SAS rmfgen and arfgen 
tools. We extracted the background from a source-free 
region. 

The object was also observed by Chandra three times 
with the net exposure time of 14 ks (seq:600114), 44 ks 
(seq:600093) and 83 ks (seq:600305), respectively. We an- 
alyzed the level 2 Chandra ACIS-S event data. We used 
the merged data of three observations without any good- 
time-interval filtering because the background count rate 
is always <6.5 % of the source in 5.0-8.0 keV. RMFs and 
ARFs were created using the CIAO mkacisrmf and mk- 
warf tools. We extracted the background from a source- 
free region. 

3. ANALYSIS & RESULTS 

3.1. Neutral, He-like and H-like Fe K lines with Suzaku 

and XMM 

Spectra taken with Suzaku were extracted from a cir- 
cular region with a radius of 2'. 3, indicated as R0 in 
Figure Q] left. The center of R0 corresponds to the peak 
position of the Suzaku 4.0-10.0 keV image. The size of 
R0 was optimized to realize the largest number of pho- 
tons avoiding other disk-origin components and anomaly 
columns in the XIS0 detector below the nuclear region. 

As shown in Figure [2] left, complex Fe K line struc- 
ture was detected from the Suzaku spectrum as large 
residuals to a power-law continuum model. The contin- 
uum was determined in 5.0-8.0 keV but excluding 6.3-7.1 
keV to avoid Fe K emission lines. We added Gaussian 
emission line models one by one. The line width was 
fixed to 0. Fitting was significantly improved by each 
step. The best-fit line center, equivalent width (EW), 
and flux, along with their errors are summarized in Ta- 
ble [1] We considered only the error of line flux in esti- 
mating the error of EW. A systematic error in flux due 
to PSF modeling is bracketed by comparing two arfs: 
one for a point source and one for a uniformly extended 
source. The difference, ~10 %, was smaller than the sta- 
tistical error. Suzaku XIS achieved an energy scale accu - 
racy better than 6 eV at 6 keV (jKovama et al.l [2007) . 
which means that the systematic error on the energy 
center is much smaller than the statistical error. The 
3 a acceptable range of the energy center of the lowest- 
energy line from our Suzaku obs ervation is 6. 29-6.43 keV, 
which includes Fe I to Fe XVI (|Housdll969r i. The most 
natural explanation of the line is the neutral fluorescent 
line. We ther e fore i nterpreted the line as Fe I, although 
iWeaver et al.1 (|2002f ) indicated that it could be emission 



from low ionization-state ions (Fe IV- Fe XV) . 

We emplo y ed the Monte Carlo procedure described in 
iBuote et al.l ([2009) to calculate a significance of these 
lines. In this process, we evaluated the reduction in the 
C-statistic (29.9 for Fe I and 11.4 for Fe XXVI) after 
adding a line with the energy center range restricted to 
6.3-6.5 keV for the Fe I line and 6.87-7.07 keV for the 
Fe XXVI line. The two lines were added individually. 
We conducted 10000 Monte Carlo simulations without 
Fe I and Fe XXVI lines with the same exposure and 
binning as the observation. The fraction of simulation 
that improved the C-statistic by 29.9 and 11.4 are and 
1.1 x 10~ 3 , respectively. Therefore the significance of the 
detection of Fe I and Fe XXVI are >99.99 % and 99.89 %, 
respectively. This is the first clear detection of the Fe I 
and Fe XXVI lines in NGC 253. 

We also extracted a spectrum from R0 by XMM as a 
complementary analysis. We found the same complex 
Fe K lines as shown in Figure [2] middle and we con- 
firmed that resultant parameters are almost consistent 
with those of Suzaku which arc summarized in Table [TJ 

3.2. Spectral Analysis with Chandra 

We utilized Chandra, characterized by a sub-arcsecond 
angular resolution, for determining the emission region of 
the Fe K lines. From a Chandra image, we see that most 
of the emission above 4 keV is concentrated in ~ 8" area 
(Figure Q] right). We divided the compact area into two 
regions, Rl and R2, as indicated in Figure Q] right to 
examine whether the emission lines are associated with 
the AGN. Rl is a circlular region of a radius of 1".2 
whose center is the peak of 4-10 k eV surface brightness. 
The b uried AGN is located in Rl (|Muller-Sanchez et all 
I2010T ) . R2 is an ellipse with the major and minor axes of 
8" and 3", respectively, but excluding Rl. First, we ex- 
tracted the spectrum from the whole (Rl + R2) region. 
The spectrum around the Fe K line complex is shown in 
Figure[2] right. We performed the same analysis as we did 
with the Suzaku data. The results are summarized in Ta- 
ble [1] Flux of Fe XXV and Fe XXVI are consistent among 
three observatories to within statistical errors, which sug- 
gests that these lines mostly originate from this compact 
region. On the other hand, the flux of Fe I observed with 
Suzaku and XMM is twice larger than that with Chan- 
dra. Hence, the neutral emission may be more extended 
than the ionized emission or variable. The smaller EWs 
of the Suzaku and XMM observations are because of a 
higher continuum due to the emission from the disk. 
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Fig. 2.— Spectra with residuals in 5.0-8.0 keV of Suzaku FI (left), XMM PN (middle) from RO and Chandra ACIS from (Rl + R2). The 
regions arc defined in Figure [l] Each residual panel indicates that to the model of power-law, (ii) power-law + gaussian, (iii) power-law + 
2gaussian and (iv) power-law + Sgaussian. The power-law continuum was determined by excluding 6.3-7.1 keV to avoid the contamination 
from Fc K emission lines. 



TABLE 1 

Energy and flux of complex Fe K lines. 



Region 




neutral Fc 


Fe XXV 


Fc XXVI 


Suzaku 




Energy (keV) * 


6.36±0.03 


6.67±0.02 


6.96_ 08 


RO 


Flux f 


3.8±i.o 


8-4±i:3 


9 O + 0.9 
"-1.2 




EW (eV) 




220^29 




XMM 




Energy (keV) * 


6.49±0.02 


6.72±0.01 


6.99±0.03 


RO 


Flux t 


, o+l.O 

J - s -o.s 


8.0±0.9 


„ ,+0.8 
■^'^-O.S 




EW (eV) 


77+20 
11 -16 


189±21 


52 + J 2 ° 2 


Chandra 




Energy (keV) * 


6.40 (fix) 


6 66 +0 - 05 

— 0.06 


6.97 (fix) 


Rl 


Flux f 


<0.9 


2.1±1.2 


<1.1 




EW (eV) 


<33 


253±145 


<91 




Energy (keV) * 


6 - 45 -0.09 


6 - 64 -o.oi 


6.97 (fix) 


R2 


Flux * 


1.2±1.1 


4 8 +1 4 
l.i 


<0.8 




EW (eV) 


86±79 


' -141 


<59 




Energy (keV) * 


6 43+ 20 


6 64 +0 - 02 

u - u ^-0.03 


6 95+ 13 


R1+R2 


Flux t 


1 7+ 11 

L - ' -1.3 


7 1 + 15 

' ■ 1 -0.9 


1.7±1.2s 




EW (eV) 


7r + 49 
'°-57 


427±°° 


82±58 




Energy (keV) * 


6 46+ 020 


6.65±0.04 


6.97 (fix) 


PI 


Flux * 


o.8i° Q ;I 


1 9+ 10 


<0.9 




EW (eV) 


176^ 


662_ 209 


<73 



* Energy center at the rest frame. 

T Photon flux in the unit of 10 -6 photons s — 1 cm -2 . 



We also extracted the spectra from Rl and R2 to eval- 
uate the contribution from the buried AGN. As is shown 
in Table [TJ Fe I and Fe XXV lines are observed in the R2 
region, while only an upper limit is obtained for Fe XXVI. 
The contribution to Fe I and Fe XXV is even larger in R2 
than in Rl. Thus the AGN is not the dominant source 
of the two lines and these emission lines are distributed 
over the (R1+R2) regions. 

3.3. Imaging analysis with Chandra 

As described above, the observed Fe K complex is un- 
likely to be explained by the AGN alone. In this section, 
we conducted imaging analysis to demonstrate the spa- 
tial distribution of this structure. Figure [3] shows the 
Chandra narrow band line flux image for the Fe XXV 



(6.6-6.8 keV) line. Contribution of the continuum was 
subtracted from the images in the following way. We 
simply assumed that the continuum above 5 keV is ex- 
pressed by a single power-law with the same photon in- 
dex (2.0) in the whole region. The contribution from 
the continuum in 6.6-6.8 keV is calculated pixel by pixel 
using the photon index and the counts in 5.0-6.0 keV, 
where no line feature exists. The assumption of a uni- 
form photon index is rather crude, because this can be 
applied for only the same spectrum properties. We also 
confirmed that the following results did not change with 
different assumption of photon indices (1.0 and 4.0). In 
Rl, a heavily absorbed non thermal component from the 
buried AGN instead of a high temperature plasma domi- 
nates the hard band above 3 keV. Despite poor statistics, 
the map suggests that the line emission is not concen- 
trated in Rl, i.e., the location of the buried AGN. In 
contrast, two bright regions of Fe XXV line are found in 
PI and P3 defined in Figure |31 

We also extracted the line flux image for the neutral 
Fe K line. However, the statistics are too poor to do any 
analysis. Hence we focus on the Fe XXV. To confirm the 
strong line emission in PI and P3, we extracted spectra 
from four regions as indicated in Figure [3] denoted as PI, 
P2, P3 and P4. PI and P4 are circles with a radius of 
1".5 while P2 and P3 are rectangles of 1".5 x 3". PI 
and P3 include the peak position while P2 and P4 are 
selected from darker areas for comparisons. PI shows a 
strong line feature, while no line feature is found in P4. 
The resulting line fluxes in PI are shown in Tabled] ~50 
% and 30 % of the observed total Fe I and Fe XXV lines 
are emitted only from PI as shown in Figure 2] left. P2 
and P3 have the same continuum flux but the resulting 
line flux of P3 for Fe XXV is twice larger than that of 
P2. 

The coordinates of peak pixel centers in PI and P3 
are (a 20 oo, 62000) = (00 h 47 m 33 s .20, -25°17'17".2) and 
(00 h 47 m 33 s .01, -25°17'19".7), respectively. These are 
coincident with the peak pos itions of 1.3 mm con tin- 
uum in the radio observation (Sak amoto et al.ll2011[) to 
within the position resolutions. The 99 % uncertainty 
of absolute X-ray position by Chandra is 0".8, while the 
beam size of the radio observation is l".lxl".l. Thus, 
the Fe XXV line shows a clear association with the two 
peaks of molecular clouds. PI is the most luminous in 
mid infrared and is suggested to be a super star cluster 
(|Keto et al.lll999D . P3 is the brightest in the 1.3 mm 
radio continuum and is featured by the water maser as- 
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sociated with star formation (e.g. lHenkel et aLl l2004). 

3.4. Wide-band Spectrum fitting of the PI region 

To extract physical parameters we fitted the wide- 
band spectra (0.8-10 keV) of PI, where the highest line 
flux was confirmed. We adopted two-temperature ther- 
mal plasmas in collisionally ionization equilibrium plus a 
Gaussian whose centroid is 6.4 keV as the neutral Fe K 
line. One plasma represents highly-ionized Si and S emis- 
sion lines and the other makes up for continuum above 3 
keV and the Fe XXV emission line. The spectrum with 
the best-fit model is plotted in Figure |U The abundance 
is fixed at solar values. The spectrum is well represented 
with the C statistic/do./ of 132/95. The temperature 
of the higher-temperature plasma is 2.5^Qg keV. This 
plasma is responsible for the Fe XXV and Fe XXVI lines. 
The unabsorbed luminosity, thermal energy and total 
iron mass of the higher-temperature plasma are calcu- 
lated to be (8.6±1.4) x 10 38 erg s" 1 , 5.4±^;| x 10 52 erg and 
20±2 Mq, respectively. The unabsorbed luminosities of 
the Fe I and Fe XXV lines in PI are also calculated to be 
1.2± J x 10 37 and 2.7±i;§ xlO 37 erg s" 1 . Note that PI 
has a radius of 24 pc and the spherical volume of 5.8 x 10 4 
pc 3 . 
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Fig. 3. — Narrow band line images by Chandra for the Fe XXV 
line (6.6-6.8 keV) after subtracting a continuum contribution. The 
color is in units of cts (141 ks) (1 pixel) -1 . The vignetting 
is not corrected for and the scale is linear. Green and magenta 
regions indicate regions from which spectra were extracted in 
this paper. Three X - ray s ources (X-l, X-2, and X-3) found by 
Muller-Sanchcz et al. (2010) are indicated by cyan. 



4.1. Point sources 

CVs and active binaries are bright point sources that 
emit Fe lines with typical luminosities of L2-10 koV ~ 
10 31 - 10 32 erg s -1 and L 2 -io keV < 10 31 erg s _1 , re- 
spectively. The luminosity in PI corresponds to 10 7-8 
sources, or a number density of ~100-1000 pc -3 . We 
can also estimate the number density from line flux of 
CVs which are relatively luminous compared to active 
binaries. The typical lumi nosity of the Fe I line i s on the 
order of IO 29 - 30 erg s" 1 (jEzuka fe Ishidalll999D . which 
leads to the large number 10 7-8 sources and a number 
density of ^500-5000 pc -3 . The number density of CVs 
and active binaries are both 10~ 5 — 10~ 4 pc~ 3 in the 
Galaxy fe.g. iRogel et al.1l2008t iFavata et al.lll995f ). Even 
at Galactic center, the expected number density of CVs 
is on the order of 10~ 2 — 10~ 4 pc -3 , assuming that their 
number density ca n be scaled from the stellar density 
(M uno et akl 120031 ). Therefore the contribution of these 
point sources seems low. HMXBs are also one of candi- 
dates with typic al luminosity of 10 33 ~ 34 erg s _1 for Fe I 
fluorescent line (|Torreion et al.l feOlO). This corresponds 
to 10 3 ~ 4 HMXBs or lO^-lO" 1 pc" 3 in PI. Considering 
the number of HMXBs in the Galaxy (^100), contribu- 
tion of HMXBs also seems low, unless extreme number of 
HMXBs and any other point sources such as black hole 
candidates are produced by star forming ac tivity, or pop- 
ulatio n of undetected obscured HMXBs ([Walter et al.1 
2006) is much larger. 

4.2. SNRs 

Supernova remnants (SNRs) likely make the biggest 
contribution as a diffuse plasma. The luminosity 
and thermal energy of young or middle-aged SNRs 
which contribute to highly-ionized Fe lines in LMC 
are 1Q 36 ~ 37 erg s" 1 and lfj 51-52 erg, respectively 
( Hug hes" et al.l [1998T ) . Both the luminosity and thermal 
energy of the higher-temperature plasma of NGC253 cor- 
respond to about 100-1000 and 10-100 young or middle- 
aged SNRs, respectively. We also estimated the number 
of SNRs using the observed Fe mass (20±2 M ) and 
Fe mass ejection pe r one Type II SN (8.4 xlO -2 M Q ; 
llwamoto et al.l fl999f ). The resultant number of SNRs 
is about 240±20 which is consistent with the estimated 
value by the luminosity. This suggests that the very 
dense massive star cluster in PI exists, which is con- 
sistent with the results of the mid infrared observations. 



4. ORIGIN OF THE EMISSION LINES 

In this section, we discuss the origin of the neutral Fe I, 
ionized Fe XXV and Fe XXVI lines. This line structure 
likely extends to the R2 region, i.e., the central 8" cor- 
responding to ^130 pc in the nuclear region of NGC253. 
However, we focus on the very compact PI region with 
a radius of 24 pc, which shows the most luminous line 
structure. We calculate the contribution from (i) point 
sources such as cataclysmic variables (CVs), active bina- 
ries and high-mass X-ray binaries (HMXBs), (ii) SNRs 
and (iii) cold molecular clouds. Candidate (i) would pro- 
duce all three lines, while (ii) and (iii) would only be 
responsible for the ionized lines and the neutral lines, 
respectively. 



4.3. Molecular clouds 

Molecular clouds that are irradiated by surrounding X- 
ray sources would emit a fluorescent Fe I line. We assume 
that the incident X-ray is photoelectrically absorbed by 
the molecular gas that subtends an angle A£l seen from 
the continuum source. The line flux of the Fe I K line 
can be estimated as 

h.i koV = 6 J dE J ds n Fc a Fc (E)F(E), 

(1) 

where e, 71f , cr Fc (E) and F{E) indicate the fluores- 
cence yield of Fe I K line (e = 0.34), the density of the 
molecular gas, the photoelectric cross section, and the 
X-ray source intensity, respectively. Assuming a ther- 
mal bremsstrahlung spectrum for the continuum source, 
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Fig. 4. — Chandra ACIS-S spectra of PI. Left: spectrum in 5.0-8.0 keV fitted with power-law + 2gaussian corresponding to Fe I and 
Fe XXV lines (red). Right: wide-band spectrum in 0.8-10 keV. Magenta, red and green lines indicate thin thermal plasmas and the neutral 
Fe K line. 



the X-ray source intensity, F(E) oc exp(— E/IcbT), 
can be determined in the 6.0-7.5 keV band. The 
Fe column density is written as J n^eds = Z-p N{r, 
where Zp e and Nh represent the elemental abundance 
of iron to hydrogen and the column density of the re- 
quired molecular gas, respectively. The photoelectric 
cross section of Fe above 7.1 keV is approximated as 
a v r (E) = 3.8 x 10-^(^/ 7.1 keV)" 2 58 cm 2 as described 
in lNobukawa et ail (|2008f ). Substituting these values, we 
obtain 

N H = (1.0±g;i) x 10 24 Ql)" 1 cm- 2 . (2) 

The solid angle ft — 2ir corresponds to the case that 
the X-ray sources are distributed entirely behind the 
molecular gas. This column density was consistent with 
the molecular clouds o bserved in the radio observation 
(| Sakamoto et al.l 120111 ) to within the statistical errors. 
Thus, the observed Fe I line flux is explained by the emis- 
sion from the NGC253 nucleus and the dense molecular 
clouds in the nuclear region. 



5. CONCLUSIONS 

Suzaku and XMM observations revealed complex Fe 
K line structure in the nuclear region of NGC 253 and 
resolved it into Fe I, Fe XXV and Fe XXVI lines. Chandra 
demonstrated that the emission is concentrated mostly 
in the ~60 arcsec 2 nuclear region, but is not localized at 
the position of a buried AGN. The Fe XXV line shows 
an association with molecular clouds. The complex Fe 
K line can be explained by a combination of 10-1000 
SNRs and a fluorescent line from molecular clouds (Nu ~ 
10 24 cm" 2 ). 
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